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We present high angular resolution [O'.'Q) observations made with the Very 



O 

! Large Array (VLA) of the radio continuum at 8.3 and 43 GHz as well as H92a 

and H53a radio recombination lines (RRLs) from the nearby (~3 Mpc) starburst 
^ I galaxy M82. In the continuum we identify 58 sources at 8.3 GHz from which 19 

^ ! have no counterparts in catalogs published at other frequencies. At 43 GHz we 

identify 18 sources, unresolved at 0''6 resolution, from which 5 were unknown pre- 
viously. The spatial distribution of the H92a line is inhomogeneous; we identify 
27 features, about half of them are associated with continuum emission sources. 
Their sizes are typically in the range 2 to 10 pc. Although observed with poorer 
signal to noise ratio, the H53a line is detected. The line and continuum emission 
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are modeled using a collection of HII regions at different distances from the nu- 
cleus. The observations can be interpreted assuming a single- density component 
but equally well with two components, if constraints originating from previous 
high-resolution continuum observations arc used. The high-density component 
has a density of ~ 4 x 10^ cm~^. However, the bulk of the ionization is in regions 
with densities which are typically a factor 10 lower. 

The gas kinematics, using the H92q; line, confirms the presence of steep ve- 
locity gradient (26 km s~^ arcsec"^) in the nuclear region as previously reported, 
in particular from observations of the [Ne II] fine at 12 iim. This gradient has 
about the same amphtude on both sides of the nucleus. As this steep gradient is 
observed not only on the major axis but also at large distances along a band of 
PA of ~ 150°, the interpretation in terms of x2 orbits elongated along the minor 
axis of the bar, which would be observed at an angle close to the inclination of the 
main disk, seems inadequate. The observed kinematics cannot be modeled using 
a simple model that consists of a set of circular orbits observed at different tilt 
angles. Ad-hoc radial motions must be introduced to reproduce the pattern of the 
velocity field. Different famihes of orbits are indicated as we detect a signature 
in the kinematics at the transition between the two plateaus observed in the NIR 
light distribution. These H92q; data also reveal the base of the outflow where 
the injection towards the halo on the Northern side occurs. The outflow has a 
major effect on the observed kinematics, present even in the disk at distances 
close to the nucleus. The kinematical pattern suggests a connection between the 
gas flowing in the plane of M82 towards the center; this behavior most likely 
originates due to the presence of a bar and the outflow out of the plane. 



1. INTRODUCTION 

M82 is an excellent candidate to investigate the physical conditions of a starburst galaxy 
because it is one of the nearest (~3 Mpc) and brightest objects of this class. The bulk of 
the ionized gas in this galaxy is located in regions that are surrounded by large amounts of 
dust: in M82 the extinction in visual magnitudes. Ay, ranges from a few to about 15 mag. 
Radio-wavelength observations, which are not affected by dust obscuration, can play a key 
role in the determination of the physical properties of the ionized gas in starburst galaxies. 

M82 has been observed in the radio continuum over a wide range of frequencies (Kron- 
berg, Biermann & Schwab 1985; McDonald et al. 2002) and is known to host a population 
of compact sources as observed at angular resolutions < 0''2 (McDonald et al. 2002 and 
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references therein). The total number of compact sources decreases at 20 cm as compared to 
2 cm; this result has been interpreted as due to free- free absorption by ionized gas in compact 
HII regions. Of the 46 compact sources identified by McDonald et al. (2002), ~ 35% were 
classified as HII regions based on their continuum spectra. 

The distribution and kinematics of the ionized gas in the central kiloparsec of M82 have 
been previously studied using radio recombination lines (RRLs). The first detections of RRLs 
from M82 (Seaquist & Bell 1977; Chaisson & Rodriguez 1977) were a major achievement 
that led to further investigations of extragalactic RRLs. VLA observations of RRLs up to 
8 GHz by Anantharamaiah & Goss (1990) toward NGC 253 motivated further intcrfcrometric 
observations of RRLs at the same frequency (Anantharamaiah et al. 1993; Zhao et al. 1996; 
Mohan et al. 2002) and higher frequencies (Zhao et al. 2000; Anantharamaiah et al. 2000) 
toward starburst galaxies. Using the total integrated line emission, global estimates were 
made for the properties of the ionized gas in M82 (Seaquist, Bell & Bignell, 1985). Using 
the Westerbork Synthesis Radio Telescope, Roclfscma (1987) obtained the velocity field at 
moderate angular resolution (~ 13") using the H166tt RRL (at 1.4 GHz). These observations 
show the rotation of the ionized gas in the central 600 pc, with solid body rotation within a 
radius of ~ 170 pc or 8". The velocity field has also been obtained using the [Ne II] line in the 
mid IR with 2" angular resolution (Achtermann and Lacy 1995). The velocity fields obtained 
from the HlGGo; and the [Ne II] line observations are not consistent, specially in the SW half. 
However, the different angular resolutions achieved for each line prevent a direct comparison. 
Higher angular resolution observations of RRLs were necessary to understand the kinematics 
of the ionized gas. Seaquist et al. (1996) observed the H41q; line with angular resolution of 
4" and, in addition to the normal rotation, they showed the presence of kinematical features 
with velocity deviations up to 150 km s^^. From observations of the ^^CO, ^'^CO and ^^CO 
lines, Weiss et al. (1999) and Matsushita et al. (2000) reported evidence for an expanding 
supershell on the SW side of the nucleus of M82. 

The high level of star formation activity at the center of M82 could have been triggered 
due to the close interaction with the neighboring galaxy M81 and the presence of a bar that 
would drive the gas inwards to feed the starburst. From observations of the morphology 
and the kinematics of the different constituents in the inner part of M82, the presence of xl 
and x2 orbits has been suggested to indicate the existence of a bar (Achtermann and Lacy 
1995; Wills et al. 2000; Greve et al. 2002). In this scenario, the ionized gas is mainly found 
along the x2 orbits i.e. highly confined near the center. X-ray observations also suggests the 
presence of a low luminosity AGN in M82 (Matsumoto & Tsuru 1999). 

The interpretation of the RRLs is not straightforward since the fine emission mecha- 
nism could involve three different contributions: spontaneous as well as internal and external 
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stimulated emission. The radio continuum emission in the starburst regions has contribu- 
tions from free-free (thermal) and synchrotron (non-thermal) radiation, which also leads to 
complexity in the interpretation of the observations. Using the radio continuum and RRL 
observations, Anantharamaiah ct al. (2000) developed a model that consists of a collec- 
tion of HII regions in order to determine the physical properties of the ionized gas in Arp 
220. These authors have been able to reproduce the observations in the frequency range 
0.15 — 113 GHz and the simultaneous existence of both low density (~10^ cm~^) extended 
(~5 pc) HII regions and high-density (~10^ cm~^) ultra-compact (~0.1 pc) HII regions was 
deduced. 

In this paper, using the Vey Large Array (VLA) of the National Radio Astronomy 
Observatory (NRAO), we present observations toward M82 of the radio continuum at 8.3 
and 43 GHz and the H92q; (3.6 cm) and H53q; (7 mm) RRLs. With an angular resolution of 
0''6 (~9 pc), we can obtain detailed information of the spatial distribution and kinematics 
of the ionized gas within the central starburst region. In particular we search for continuum 
emission sources associated with the H92q; and H53q; line emitting regions. This paper 
is organized as follows: the observations are discussed in § 2; the results are presented in 
§ 3; a model, based on the observations of radio continuum and RRLs, that consists of a 
conglomerate of HII regions is presented in § 4; a discussion of the radio continuum emission 
at 8.3 and 43 GHz, the implications of the kinematics and the results of the proposed model 
are discussed in § 5; and finally, the conclusions are summarized in § 6. 



2. OBSERVATIONS. 
2.1. 8.3 GHz data. 

2.1.1. Observations and calibration 

The observations of the H92a line {urest =8309.3832 MHz) were conducted in the C 
(April 22, 1996), CnB (Feb 13, 2000) and B (May 05 and 11, 2001) VLA configurations. The 
maximum angular resolution achieved is ~ 0''6. The observations made in the C array have 
been acquired from the VLA archive database. The 31 spectral channels mode was used 
with a total bandwidth of 25 MHz, corresponding to a velocity coverage of ~850 km s~^, 
centered at a heliocentric velocity of 200 km s~^. The flux density scale was determined by 
observing 3C286 (5.3 Jy at 3.6 cm). The phase cahbrator was 1044-1-719, with a flux density 
of 1.5 Jy. The bandpass calibration was made using the calibrator 3C48 with a flux density of 
3.2 Jy. Bandpass calibration is critical in these observations because the line-to-continuum 
ratios are < 1%. The data were Hanning-smoothed offline to improve the signal-to-noise 
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ratio (S/N) and minimize the Gibbs effect. Tfie effective velocity resolution is 56 km s~^ 
(806 kHz). Each of the databases were self-calibrated in phase using the continuum channel; 
the solutions were then applied to the spectral line data. The data sets taken at different 
epochs were combined into a single data set. The self-calibration process was further used 
on this data set to correct for small phase offsets between the individual observations. The 
AIPS task UVLSF was used to estimate the continuum level by fitting a linear function 
through the spectral channels free of line emission, and then subtracted from each visibility 
record. The observational parameters for the H92q; hne are listed in Table 1. 



2.1.2. Imaging. 

The use of a natural weighting scheme for the data in the C configuration relative to 
those in the B configuration produces a synthesized beam with prominent wings (at the 
level of 2%). In order to correct this problem, we use a technique that consists of re- 
weighting the data in the u,v plane so as to produce a gaussian beam in the image. After 
the re-weighting process, if necessary, the image is deconvolved using the CLEAN algorithm. 
Because the u,v plane is properly sampled in the inner region, the deconvolution is not 
required when producing images at an angular resolution of 2" FWHM. For the images 
made at angular resolutions of O'.'G and 0''9 (section 3), a deconvolution is required. In these 
cases the 2" angular resolution image is used as an input for regularization in the process 
of deconvolution. The quality of the final images is limited by errors due to the imperfect 
continuum subtraction. These errors appear as weak fluctuations over a scale of one arcmin. 
These fluctuations have a typical correlation length of about 4 MHz along the spectral axis. 
This behavior results in systematic errors which hmit the quality of the "baselines" (~ 2%) 
in the spectra produced from different regions of M82 for our data analysis. 



2.2. 43 GHz data 

2.2.1. Observations and calibration 

The observations of the H53q; line {urest =42951.9714 MHz) were carried out in the 
C array of the VLA on April 13, 2000. The VLA correlator is limited to a bandwidth of 
50 MHz (~350 km s^^ at 7 mm). At 43 GHz three adjacent spectral windows, each of 
them with 15 channels, arc required to cover the velocity range of the RRL (~600 km s^^). 
The windows were centered at 42185.1, 42214.9 and 42235.1 MHz. The amplitude, phase 
and bandpass calibrators used were 3G286 (1.47 Jy), 0954+658 (0.87 Jy) and 1226+023 
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(19.4 Jy), respectively. The flux density calibrator was observed only for the central LO 
window (42214.9 MHz). The flux densities of the phase and bandpass calibrators in the 
adjacent LO windows (42185.1 and 42235.1 MHz) were assumed to have the same values 
as those determined from the central LO window. The phase calibrator was observed at 
time intervals of 10 min. The on-source integration time was ~2 hrs for each frequency 
window. The bandpass response of the instrument is different for each frequency window; 
the observation of both the bandpass and the phase calibrator must be interleaved between 
the different overlapping frequency windows required to observe the complete line. This 
method removes the offsets (< 5%) between the frequency windows. In order to correct 
for phase decorrelation that may be present for time intervals <8 min, a phase calibration 
was performed initially followed by a second calibration step applied to both phase and 
amplitude. Phase correction is more important at these frequencies since the troposphere 
introduces phase offsets that may affect the coherence of the data. The line-to-continuum 
ratio for this line is 10 — 20% in the brightest region of M82. Table 1 shows the details of 
the observations for the H53q; line. 



2.2.2. Imaging. 

In order to produce the H53a line image, the contribution of the continuum emission 
must be removed and the three frequency windows combined. The method uses two a 'priori 
constraints on the source model that are not independent. The flrst constraint is set by 
assuming that the emission in the lowest spatial frequencies is mostly continuum emission, 
allowing us to clip the data for the shortest spacings in the u,v plane. The second constraint 
is based on the fact that the projected velocity of the gas in M82 varies across the image 
due to the rotation of the galaxy. Hence, if we include a priori information on the rotational 
properties of M82, for each position in the image, we can maximize the number of line-free 
channels in each LO setting to reduce the uncertainties when removing this contribution. As 
a result this second method involves the following steps for each LO setting: (1) filter out 
part of the emission present at the lowest spatial frequencies, (2) produce an undeconvolved 
3-dimensional image, (3) determine the level of remaining continuum for each position, (4) 
determine the offsets by comparing these levels from one LO setting to the other for each 
position, (5) reflne the determination of this remaining continuum by including the line- 
free channels of the adjacent LO settings for each position, (6) subtract the determined 
continuum contribution and (7) average the three spectral windows to make the single 3- 
dimensional line image. It is not necessary to deconvolve the results after step 6 because the 
signal-to-noise ratio for the line emission is between 2 and 5. 
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In order to produce the 43 GHz continuum image we carried out the following steps: (1) 
produce a raw 3-dimensional image for each LO setting, (2) subtract the line contribution as 
derived above in each LO setting, (3) determine the mean continuum level for each position 
in each LO setting, (4) average the three 2-dimensional images obtained in step 3 and also 
average the three synthesized beams and (5) finally deconvolve this averaged image. 

The undeconvolved images were generated using the AIPS IMAGR program applying 
the suitable u,v tapers to obtain images at different angular resolutions. They were further 
processed in the GIPSY environment to perform all the subsequent steps. The data were 
Hanning-smoothed offline and a final spectral resolution of 44 km s~^ was achieved. 

3. RESULTS 

Globally, from the 8.3 GHz data, we obtain a total integrated continuum flux density of 
2.6±0.1 Jy and a H92« line flux density of 2.9±0.1 Jy km s~^. At 43 GHz the total continuum 
flux density is in this case 0.82±0.16 Jy and the H53a line flux density is 3.0±0.2 Jy km s^^. 
Images at three different resolutions are produced to analyze the data. The high angular 
resolution (0''6) images are suitable to investigate the compact bright features. Using the 
intermediate resolution images (0''9), it is possible to study the weak extended features with 
sufficient angular resolution. The low angular resolution images (2") provide information on 
the overall distribution of the ionized gas. 

3.1. Radio continuum at 8.3 and 43 GHz 

A correction for the primary beam was applied to the images at 8.3 GHz and 43 GHz. 
In the images at 8.3 GHz this correction does not exceed ~2%, while for the 43 GHz images 
this correction is a factor of 1.2 at the extreme edges of M82. The radio continuum images 
at 8.3 and 43 GHz are shown in Fig. la and Fig. lb, respectively; the high angular resolution 
images (contours) are shown overlaid on their corresponding low angular resolution images 
(gray scale) at each frequency. 

A number of previous high angular resolution (< 0.1") studies carried out over the 
range 408 MHz to 15 GHz have revealed the existence of a population of compact sources 
classified as HII regions or SNRs (Kronberg et al. 1985; Huang et al. 1994; Muxlow et al. 
1994; Wills et al. 1997; Allen & Kronberg 1998; McDonald et al. 2002). In order to study the 
small scale continuum features at 8.3 GHz, the underlying extended emission was filtered 
out by restricting the lower uv-baseline range (to u,v distances larger than 50 kA) so that no 
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structures larger than 4" are present. At 43 GHz an iterative procedure is used to determine 
and remove the spatiaUy extended continuum emission. Figures 2a and 2b show the resulting 
spatial distribution of these features at 8.3 and 43 GHz, respectively. In these images we 
are able to identify ~ 85% of the previously known compact sources. The remaining 15% 
of compact sources were not detected because their peak flux density is below the level of 
the local diffuse emission. On the other hand, these images reveal 21 new features (19 of 
these sources are detected at both frequencies and two sources only at 43 GHz) not detected 
in previous studies. Several of these sources are relatively more extended (> 0''5) and they 
could not be identified in the previous high angular resolution images. Including our new 
detections, a total of 60 sources are detected toward M82, 58 features at 8.3 GHz and 18 at 
43 GHz. 

In Table 2 we list the parameters of these 60 compact features: 

• The first column hsts the galactic names for each source. For each feature we measure 
its spatial position, its deconvolved angular size and flux density. The measurements 
at 8.3 GHz are summarized in columns 2, 3 and 4 and those at 43 GHz in columns 5, 6 
and 7, respectively. For the 21 newly detected features a super-script "a" is used. The 
uncertainties in flux densities correspond to a 1 arms noise in the least-squares fitting 
procedure. For features not detected, we provide an upper limit of 3 arms for the flux 
densities. 

• From the measured flux densities (columns 4 and 7) we derive the spectral index 
(column 8) between 8.3 and 43 GHz. This value can be compared with the spectral 
index previously reported in the literature (column 9) using different frequency ranges 
(e.g. 0.408 to 15 GHz). 

• In column 10 we list identifications for the sources which have radio recombination line 
emission. The associations arc made with the H92a line features observed at 0"9 (see 
section 3.2). At an angular resolution of 0!'6, only four H92q; compact line features are 
detected with signal-to-noise ratio > 3 (see section 3.2 and Fig 3). 

• In the last column, the classification is given (HII regions or SNR). The sources previ- 
ously identified as HII regions by McDonald ct al. (2002) were classified based on their 
spectral indexes, a, which in these cases arc inverted {S^ oc z/", a > 0) between 5 and 
15 GHz (column 9). We confirm the classification for all the HII regions previously 
reported, as they have a spectral index that is consistent with optically thin free-free 
emission between 8.3 and 43 GHz (column 8). We assign tentatively an HII type to 
four more features, because between 8.3 and 43 GHz, 43.21-1-61.3 and 45.63-1-66.9 have 
flat spectra and 44.17-1-64.4 and 44.43-1-62.5 have inverted spectra. For the sources 
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previously identified as SNRs, our results show non-thermal spectra between 8.3 and 
43 GHz (sec column 8) as expected from the spectral index computed between 5 and 
15 GHz. In addition to those features previously classified as SNRs, based on the upper 
limits obtained from the spectral index between 8.3 and 43 GHz, we can tentatively 
consider 39.92+55.9, 40.49+57.4, 43.00+59.0, 44.11+64.3 and 46.74+69.7 as SNRs. It 
is noted that three SNRs, 41.29+59.7, 43.00+59.0 and 45.79+65.2, exhibit H92a RRL 
emission; this issue will be discussed in § 5. 

3.2. Radio Recombination Lines H92q; and H53q; 

Figure 3 shows the high-resolution (0''6) H92a line images in the region where H92a 
RRL emission is detected; in these images four sources are clearly identified. The center of 
this region is located ~4" W of the 2.2 /xm peak (Lester et al. 1990) and corresponds to the 
brightest area observed in the H92q; fine in M82 on a larger scale. All four fine emitting 
sources are spatially associated with compact continuum sources previously identified by 
McDonald et al. (2002). No other fine features are detected outside this region at this 
resolution. For Figure 4 we show the velocity-channels at an intermediate angular resolution 
of 0''9 (top) and an image of the total integrated H92« line emission (bottom). In this case 
27 line features can be identified and those associated with continuum sources are listed 
in column (10) of Table 2. In these images, the position of the HII regions (crosses) and 
SNRs (stars) are shown as reported by McDonald et al. (2002). In this comparison the 
SNRs are mainly located at the periphery of the line emitting regions while the HII regions 
are clearly associated with the RRL emission. Because these HII regions are the compact 
continuum sources which have inverted spectra between 5 and 15 GHz (section 3), we have 
direct evidence that this population of compact objects is associated with the more extended 
HII regions which produce the H92q; recombination line emission. Figure 5 (top) shows the 
H92q; line velocity-channels at a low angular resolution of 2" (top) for the overall distribution 
of the integrated H92q; line emitting regions (bottom) . 

Figures 6a and 6b show the low-resolution images of the integrated H92q; and H53q; 
line emission (contours) overlaid on the low-resolution images of radio continuum at 8.3 and 
43 GHz (gray scale), respectively. To produce these images for the integrated line emission 
only the regions with line emission > 3aj.ms were selected. There is a clear correspondence 
between the spatial distribution of the line and the continuum emissions. In both cases the 
distribution of the total line emission is characterized by a pair of concentrations on each side 
of the center defined by the 2.2 /xm peak (Lester et al. 1990). However, the sources located 
on the E side are fainter compared with their counterpart on the W side. This asymmetry 
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has also been observed in [Ne II] and for several molecules (Achtermann and Lacy 1995). 
Figure 6c shows the velocity integrated line emission of the H53q; line (contours) overlaid 
on the corresponding image of the H92q; line (gray scale). A good correspondence between 
the spatial distributions for the two lines is observed. Both lines cover an angular size of 
about 30" (480 pc) along the major axis and the main peak is ~ 4" W of the 2.2 fim peak. 
However, we note that the H92a line emission may be more extended than the H53q; line 
emission (~ 4") on the extreme NE side of M82. More sensitive H53q; line observations 
would be required to verify this possible difference. In order to characterize the physical 
properties of the ionized gas in M82 (section 4.2), we define five regions ( El, E2, C, Wl 
and W2) as indicated by the rectangular boxes in Figure 6c. Regions El, Wl and W2 have 
been labeled according to peaks observed in [Ne II] emission (Achtermann and Lacy 1995). 
The comparison of the distribution of the H92q; line emission with the distribution of other 
tracers of the ionized gas is discussed in § 5.2. 

Figure 7 shows the spatially integrated H92q; and H53a spectra for four distinct regions 
( El, E2, Wl and W2); the spectrum for region "C" is not shown since H53q; line emission 
is not detected above 3arms- The area for the integration of line and continuum emission 
are defined by the first contour (57 Jy beam~^ m s~^) of the H53a distribution. The H92« 
and H53q; line parameters derived from the fits of gaussian profiles (Figure 7) are given in 
Table 3. Columns (2), (3) and (4) of Table 3 list the parameters obtained for the RRL H92q; 
and column (5) the continuum flux density at 8.3 GHz. Columns (6), (7) and (8) fist the 
parameters obtained for the RRL H53q;, column (9) the continuum flux density at 43 GHz, 
column (10) the maximum filling factor for the high-density component explained in § 4 
and column (11) the continuum fiux density at 5 GHz (personal communication, A. Pedlar). 
These measurements form the basis for the models discussed in § 4. 

The global properties for the five regions identified in Figure 6c are given in Table 
4, the total H92q; line emission integrated over each rectangular box (column 2) and the 
corresponding line (H92q;) to continuum (at 8.3 GHz) ratio (column 3). Inside these four 
regions a number of sub-regions have been identified in the H92« line image (see Table 4). For 
each of these sub-regions (labeled in column 4), Table 4 lists the spatial coordinates (column 
5), the velocity- integrated line emission (column 6), the FWHM angular size deconvolved to 
correct for the 0''9 beam broadening (column 7) , the peak flux density (column 8) , centroid 
velocity (column 9) and the spectral width (column 10). The angular size of each sub- 
region is based on fits of 2D gaussians in the appropriate spectral channels; the spectral line 
parameters were obtained fitting a gaussian profile along the spectral axis for the spatially 
integrated emission in each of these sub-regions. For each region the sum of the line fiux 
densities of the sub-regions is also given and can be compared with the total line fiux density 
in column 2, determined from the low angular resolution image. For all regions, essentially 
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all the line flux density arises from these sub-regions. Thus, any diffuse emission is < 10%. 

The velocity field determined using the H92q; line image at 2" angular resolution is 
shown in Figure 8a. In this image, the major and minor axis are indicated for a disk with a 
diameter of 40" as observed at an inclination of 81°, and position angle (PA) of 68° relative 
to the observer. The intersection of these two axes is set at the position of the 2.2 //m peak 
(Dietz et al. 1986), which also coincides with the center of the hole in Fig. 5b. The position 
of the center of the ring fitted by Achtermann and Lacy (1995), assumed to be the kinematic 
center, is indicated by a cross. The important features in this velocity field are the following: 
a) The major axis corresponds to the line of nodes which gives the maximum velocity gradient 
for the outer regions, b) in the central region there is a significant tilt of the iso- velocities 
relative to the PA of the minor axis, c) this velocity field is dominated by an axisymmetric 
pattern and the position of the center of symmetry appears to have a small offset (~ 1''5) 
relative to the position of the 2.2 /im peak, d) Some areas show significant deviations relative 
to this general pattern, in particular in the region near 9'^51™43f5, +69°55'04", where the 
velocities are ~ 40 km s~^ larger than expected. The kinematics of the ionized gas are 
discussed in § 5.3. The spatial distribution of the velocity dispersion is shown in Figure 8b. 
One important feature in this image is the presence of a band ~ 1''5 wide with a PA of 
~ 150°, close to the minor axis, where the spectra are systematically narrower. This feature 
is analyzed in section 5.3. 



4. Models with multiple HII regions. 

In three starburst galaxies (Arp 220, M83, NGC 2146) a model consisting of a collection 
of HII regions has been used successfully to explain the observations of RRLs and the radio 
continuum (Anantharamaiah et al. 1993; Zhao et al. 1996; Anantharamaiah & Goss 1996; 
Anantharamaiah et al. 2000). The ionized gas distribution in the center of M82 is inhomoge- 
neous as observed in the RRLs H92q; and H53q;. Thus, a model that consists of a collection 
of HII regions seems appropriate to explain the RRL and radio continuum emission from this 
galaxy. Based on this model we can derive the physical parameters of the HII regions (e.g. 
electron temperature and volume density). The model must reproduce both the observations 
of the line and continuum emission as a function of frequency. For M82, observations at 2" 
angular resolution of RRLs are available at two different frequencies (8.3 and 43 GHz) and 
observations of radio continuum at three frequencies (5, 8.3 and 43 GHz). The 5 GHz flux 
densities are listed in column 11 of Table 3. 

The gas is ionized by young massive stars (types O and/or B) radiating large amounts 
(> 10^^ s~^) of Lyman continuum photons. The number of Lyman continuum photons 
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determines the size, I, of each HIT region depending on the electron density. Since these 
models are constructed to obtain average values for the physical properties of the ionized 
gas, we have assumed: (1) Each ionizing massive star emits the same number of Lyman 
continuum photons (N^^^c = 10^^ s~^) and (2) N^y^ is related to the electron density and the 
size of each Hll region by Niyc oc Ug/^; this relation implies that for each density value, / will 
be determined and then / is not a free parameter in the model, providing a constraint. Since 
massive stars have a relatively short lifetime on the main sequence, the derived number of 
Lyman continuum photons can be related to the star formation rate (SFR) of O and B stars 
if the absorption of Lyman continuum photons by dust inside the HIT regions is neglected. 

This model consists of Hll regions located in front of a diffuse mixture of thermal and 
non-thermal emission {Scbg)- The HIT regions radiate free- free emission and recombination 
lines. Background non-thermal emission and free-free emission that arise inside each HIT 
region may stimulate the emission of RRLs in the HIT regions. Thus, the line emission from 
each HIT region has three different contributions: (i) spontaneous, (ii) internally stimulated 
and (iii) externally stimulated emission. Each HII region is characterized by an electron tem- 
perature Te (K), electron density rie (cm~^), and linear size / (pc). The thermal continuum 
flux density, Sc-th, is given by, 

Sc-TH = fimiB^il - e-^^) mJy, (1) 

where Qhii is the solid angle subtended by each Hll region, B^y is the black-body emission 
(mJy) and Tc is the continuum optical depth, which depends on frequency as follows (Bell 
& Seaquist 1978): 

TC = 3.01 X W-'iy-X-'-Hn( ^ )eMc, (2) 

where EMc = n^l is the emission measure for each Hll region. The line flux density, 
Sl (mJy), arising from each HII region is calculated using. 



+ Scbge-^^{e-^^ -l)mJy. (3) 

The spontaneous and the internally stimulated emission are represented by the first 
term in Eq. 3 and the externally stimulated emission is represented by the second term, k is 
the Boltzmann constant, c is the speed of hght, t£ is the peak line optical depth under LTE 
conditions, tl — hnPn'Tl i^ peak fine optical depth corrected for non-LTE effects, 6„ is 
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the departure coefficient, n is the quantum number, and /3„ = 1 — {kTf./hvi)dlnhn/dn. For a 
given combination of parameters Tg and Ug, we compute the corresponding values of 6„ and 
Pn (Salem & Brocklehurst 1979) which are used to calculate the integrated line flux density 
of a single HII region using Eq. 3. The line optical depth is given by t£ — TcSl/ Sc-th where 
the line to thermal continuum ratio in LTE is defined as. 



where EM^ is the emission measure from hydrogen and EM^; is the emission measure from 
all ions [cmT^ pc). In these models, EMj^/EM^; = 1; for simplicity we do not consider the 
presence of ionized He. 

Single-density (SD) models, consisting of a number of HII regions with the same physical 
properties, have been used to reproduce the observations. Additionally, following Ananthara- 
maiah et al. (2000), two-density (TD) models that consist of both low- and high-density HII 
regions were used to reproduce the observed values (listed in Table 4) for four regions ( El, 
E2, Wl, W2); in region 'C the H53q; line emission is not detected (see Figure 6c). The 
following constraints have been imposed: (i) the predicted thermal radio continuum flux 
densities should not exceed the radio continuum flux density observed at 43 GHz, (ii) the 
area flUing factor should be < 1, (iii) the peak line emission from the model should not be 
larger than the corresponding observations, (iv) the spectral index of the background radi- 
ation has been limited to be not steeper than —1.5, (v) the continuum computed from the 
model must reproduce all continuum observations at 6 cm, 3.6 cm and 7 mm, and (vi) the 
area filling factor for the high-density component has an upper limit imposed by observations 
of compact continuum sources at high angular resolution (McDonald et al. 2002). For the 
two-density models the constraint (vi) is based on the additional assumption that the high- 
density HII regions are coincident with compact sources previously identified as HII regions. 
Using the emission measure given in Table 4 of McDonald et al. (2002), we compute the 
angular sizes and the maximum filling factor for the HII regions that form the high-density 
component. The models that are consistent with these constraints were accepted as valid 
solutions. 

The number of HII regions for the single-density model is computed by dividing the ob- 
served integrated line fiux density by the integrated line fiux density from a single HII region. 
In the case of the model with two-density components, the number of HII regions with high- 
density gas {J^hii-hd) is computed dividing the area filling factor derived above by the area 
filling factor that a single high-density HII region can occupy (following ng oc /~"^/^). Once 
we know the contribution of the high-density component to the line and continuum emission. 
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the number of HII regions in the low-density component (^hii-ld) can be estimated (see 
Table 5) . 

For both the single-density and the two-density models, the contribution from the back- 
ground emission is computed as the difference between the observed radio continuum at 
8.3 GHz and the value predicted by the models for the thermal emission at this frequency. 

The spectral index value (a), determined for this background emission, is obtained by mini- 
mizing the difference between the observed and the predicted values of the thermal continuum 
flux density at 43 GHz. 

Since we have assumed that each HII region is being excited by stars that emits ^Lyc — 
10^^ photons per second, the total rate of emission of ionizing photons, ^Lyc-tot, is, 

Nhyc-tot = NLycNnii, (5) 

where Nhii = Nhii-ld + Nhii^hd- Another way to estimate ^Lyc-tot is using (Schraml & 
Mezger 1969; Rodriguez et al. 1980), 

where we can use the flux density from thermal emission, Sc-th, at any frequency (if 
Tc << 1) and D, the distance to the galaxy. In the case of optically thin emission at 43 
GHz, both estimates of '^Lyc give similar values. However, the model is based on a number of 
assumptions. On the other hand, the number of Lyman continuum photons is proportional 
to the radio continuum at 43 GHz and does not depend on extinction estimation as occurs for 
the near-infrared case. Thus, the determination of the number of Lyman continuum photons 
from the continuum at 43 GHz is the most reliable way to estimate the Lyman continuum 
photons production. If the gas is optically thick (which may be the case for a high-density 
component) at 43 GHz, then the estimation of the total number of ionizing photons using 
Eq. 5 is a lower limit. The total SFR can be derived from the Lyman continuum luminosity 
following Anantharamaiah et al. (2000). 

4.1. Results from the models. 

The parameters Tg and were considered to lie in the ranges 5000 — 10, 000 K (Garay 
& Rodriguez 1983) and 500-10^ cm~^, respectively. The low-density component was deflned 
to consist of HII regions with electron density, Ug in the range 100—10^ cm~^, while the 
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high-density component was defined as composed of HIT regions with ng in the range 2 x 
10^ — 10^ cm~^. In Table 5 we fist the range of physical parameters estimated using the 
models for each identified complex of HII regions (see Fig. 6c) . For each identified complex 
of HII regions, Table 5 lists: the size of each HII region, number of HII regions, factor of 
area covered by each density component, number of ionizing photons, star formation rate, 
line and continuum opacity at 8.3 GHz, the departure coefficients for the H92q; line, the 
spectral index of the background emission and the total mass of ionized gas. Figure 9 shows 
the expected variation of radio continuum and the integrated RRL strength as a function of 
frequency for regions El, E2, Wl and W2. The models shown in Figure 9 correspond to the 
mean values based on the models hsted in Table 5. 

The single- density models reproduce the observed line and continuum flux densities for 
all regions considered. However, these results only provide a rough estimate for the physical 
properties of the ionized gas. As we know from previous studies (McDonald et al. 2002), 
there are compact HII regions in which higher density gas (> 10^ cm~^) could be present. 
Using the area filling factor constraint in the two-density models (see section 4). we obtain 
valid solutions for regions E2, Wl and W2. These two densities are a few 10^ cm~^ and 
10"^ cm~^ as listed in Table 5. For the region El, there is no evidence of any compact high- 
density HII regions (McDonald et al. 2002). Since for region El we do not have an area 
filling factor constraint, an estimate of the physical properties of the high-density ionized 
gas was not computed. 

Based on these models the area filling factor from all HII regions, including the cases 
where the high-density component is taken into account, is ~ 0.3. We consider that the 
number of HII regions (obtained from the models) in the compact features is an upper limit. 
As hsted in Table 5, the low-density ionized gas has an electron density of ~ 5 x 10^ cm~^. 
McDonald et al. (2002) found that the turnover in the spectrum of SNRs at 1.4 GHz requires 
foreground emission measures of ~ 10^ pc cm~^, consistent with our model. However, it is 
likely that there is also ionized gas in a diffuse component with lower electron densities. 
McKcith et al. (1995) estimated the electron density to be typically 1.5 x 10^ cm~^ in the 
plane of M82 near the center of the galaxy. In order to explain observations in the far infrared 
(made at an angular resolution of 80"), Colbert et al. (1999) presented a model which is a 
combination of HII regions with Ug ~ 0.25 x 10^ cm~^ and photo-dissociated regions (PDRs) 
with ~ 2 X 10^ cm~^. High resolution observations of RRLs of higher quantum number 
(e.g. H116q;), with an angular resolution similar to the H92q; data presented here, are 
required to constrain the models to account for the gas at lower density. 

Assuming that the relative number of high-density HII regions compared to the number 
of low-density HII regions is an age indicator of recent star formation activity, based on the 
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two-density models, we infer that the region W2 is the youngest. Another age indicator for 
the starburst in each region is the spectral index a: the steeper the spectral index the more 
evolved the star formation activity, which is true if we assume that each region has formed 
stars at a constant rate, at least during the last few 10® yr. Using this criteria E2 would then 
be the most evolved starburst. Observations of low quantum number RRLs (e.g. H43q;) 
tracing the dense gas are required to better constrain the recent star formation history. 

The number of Lyman continuum photons required to ionize all four complexes, ~ 
16 X 10^^ s^^, is a factor three lower than the value estimated by Achtcrmann and Lacy 
(1995). The sum of the H92a line flux densities from the four modeled regions represents 
only ~50% of the total line flux density; our analysis is thus restricted to the areas where 
the H53q; line is detected. Considering this fact, the Lyman continuum rate obtained from 
the models is a lower limit for the total Lyman continuum rate emission in the center of 
M82. Then, we conclude that our results are in reasonable agreement with previous results. 

M82 differs from Arp 220 in the sense that very high-density HII regions are not required 
to explain the hne and continuum emission from the center of M82. The SFR depends on 
the initial mass function and mainly on the upper limit used for the mass of the stars that 
are formed. Using the Miller-Scalo IMF and mass hmits of 1 and 100 M©, the total SFR 
derived for Arp 220 is ~ 240 yr~^ (Anantharamaiah ct al. 2000), whereas that for M82 
is ~ 100 times lower (~ 3 M0 yr~^). Another estimate of the SFR for M82 using the same 
mass limits, the same IMF and assuming that all the radio continuum emission is thermal 
gives an upper limit of 46 Mq yr~^ (Kronbcrg ct al. 1985). If the Salpeter IMF is used then 
the total SFR in M82 would be ~ 3 times larger. The present mass of gas in the central 
region of M82 (160 pc) is estimated to be ~ 2 x 10^ M©, a factor 30 lower compared to 
Arp 220. 



5. DISCUSSION 

5.1. Radio Continuum and RRLs H92q; and HSSa. 

There are two possible reasons that could explain why some continuum features were 
detected only at 8.3 GHz: (i) the features are HII regions and the higher noise level at 43 GHz 
limits the detection, (ii) the features are SNR with non-thermal spectrum (a < 0, S'jy oc z/") 
and flux density at 43 GHz < 3o",.ms- For all the features idcntiflcd as SNR, the spectral 
index value determined from our observations is in good agreement with a non-thermal 
spectra. From 17 sources idcntiflcd as Hll regions, 13 sources have spectral index vahics 
in agreement with a thermal spectrum ( 8.3 — 43 GHz). The compact sources 42.69+58.2, 
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42.56+58.0 and 42.48+58.4 are observed in a larger complex region; because this complex 
region is embedded in a more extended emission region with a similar flux density level, 
we were not able to determine the continuum flux densities for each of them individually. 
The spectral index from this complex region is not in agreement with a thermal spectra; the 
contribution of synchrotron emission from the SNR at 8.3 GHz may explain this result. Due 
to the lower angular resolution of our observations, the flux densities of some features are 
larger than previously reported measurements at 3.6 cm based on higher angular resolution 
data (Huang et al. 1994; Allen & Kronberg 1998). The larger measured flux densities arise 
from extended features that were resolved out by the higher angular resolution observations. 
There are three compact sources that have been classified as SNRs and have detectable 
RRL emission. The H92a line emission that seems to arise from SNRs may be accounted 
for if a group of HII regions is located along the same line of sight as the SNR. The RRL 
emission that arise in the foreground ionized gas could even be stimulated by the external 
emission arising from the SNR. Quite possibly the externally stimulated emission plays an 
important role in these lines of sight. For the source 41.95+57.5, the flux density at 8.3 GHz 
is 13.6+0.4 mJy using observations made in 2001. Allen & Kronberg (1998) at 8.4 GHz 
determined a value of 26.76+0.5 mJy in 1994-1995. The implied decay rate of this radio 
supernova is in good agreement with the estimates (8.8% yr~^) of Allen & Kronberg (1998). 



5.2. Compcirison with other tracers of ionized gas. 

As observed at an angular resolution of 2", the structure of M82 is very similar in the 
H92q; line (Fig. 5) and the [Ne II] line (Fig. 1 of Achtermann & Lacy 1995). In both cases the 
most prominent source is Wl (see Fig. 7), located ~ 5" E of the 2.2 //m peak. In the H92q; line 
images (Fig. 5), a quasi-circular "hole" with a FWHM deconvolved size of ~ 2" is observed. 
The center of this hole is located at a = 09''51"'43?5, 5 = +69°55'00M, which is close to the 
2.2 /xm peak determined by Dietz et al. (1986) at a = 09'*51™43f5, S = +69°55'00''7 and 
shifted from the one determined by Lester et al. (1990) at a = 09'^51™43f6, 6 = +69°55'00''l 
by l'.'5 in RA and by 0'.'3 in Dec. The position of this peak needs to be conflrmcd to 
deflnitively conclude if this hole as observed in H92q; is, in fact, centered on this peak. The 
regions El and Wl are connected by a faint ridge of emission on the N side of this hole (see 
Figure 6). Based on the [Ne II] observations, Achtermann and Lacy (1995) have suggested 
the presence of an ionized ring with a projected axial ratio of about three, a major axis of 
11" at a position angle of 70° and centered at a = 09''51'"43f4, 6 = +69°55'00'.'l. Since 
the emission of H92q; line is not affected by extinction, the lack of emission on the S side 
indicates that the ionized gas is not uniformly distributed along this ring. In agreement 
with the conclusions of Achtermann and Lacy (1995), the observed brightness of El and Wl 
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cannot be explained simply by a limb brightening effect. 

At 0''9 resolution (Fig. 4) the spatial distribution of the H92q; hne is quite patchy. The 
region Wl is clearly decomposed into two main components: the group composed of Wlb 
and Wlc is the brightest and is located at the SE side and the group composed of Wle and 
Wlf is located at the NW side. These two groups would be located on the S and the N 

side (respectively) of the major axis (PA 70°) of the ring proposed by Achtermann and Lacy 
(1995). These two groups, as observed in the Br'-j (2.17 /xm) line, are separated by a ridge 
of extinction that extends along the major axis of M82 (Larkin et al. 1994). Since the H92q; 
line is not affected by extinction, the lack of emission between these two groups on the W 
side of the 2.2 /xm peak is not due to dust extinction. On the E side, as observed in the H92q; 
hne, El appears relatively more diffuse than Wl, with the same extension as the Br^ line. 
This region is located on the S edge of a ridge of extinction that extends along the major axis 
of M82. As we do not detect significant H92q; hne emission on this ridge, the distribution 
of Br^j line emission for the region El does not seem to be significantly affected by dust 
extinction. The extinction derived from the Br7 line was used to determine the orientation 
of M82. However, the distribution of the Br'-y line emission is not severely affected by dust 
extinction. Thus, the fact that the brightest Br'-j sources appear on the S side cannot be 
used as a reliable indicator to determine the orientation of the galaxy (Larkin et al. 1994). 



5.3. Kinematics of the ionized gas. 

The kinematics in the central part of M82 have been extensively studied using observa- 
tions of molecular lines, HI in absorption and the infrared lines Br^ and [Ne II]. From these 
studies a picture emerges in which there is an inner ring of ionized gas ~ 11" (~ 200 pc) 
in diameter possibly surrounded by a ring of molecular gas. An alternate model for the 
molecular distribution with two lobes has also been proposed by Larkin et al. (1994). The 
presence of a stellar bar, about one kiloparsec in length, is suggested by the 2.2 /xm hght 
distribution (Lester et al. 1990). The ring of ionized gas could be related to x2 orbits at 
the center of the galaxy; observationally it is unclear if this ring is circular or has an oval 
distortion, since M82 is viewed highly inclined (i=80°). Achtermann and Lacy (1995) have 
shown that their data are consistent with a simple circular ring of about 150 pc diameter 
with a rotation velocity of ~ 112 km s^^ and possibly an associated outflow. They also 
present a model involving xl orbits with an almost end-on view of elliptical x2 orbits. 

Along the major axis of M82 the 2.2 yum light distribution is characterized by two 
plateaus, an inner plateau extending over ~ 20" and an outer one extending over one arcmin; 
the outer plateau has been interpreted as a possible evidence for a stellar bar (Telesco et al. 
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1991). Using images of M82 at the I, J, K and L' bands, Larkin et al. (1994) inferred that a 
dust lane is in front of the stellar population to the W of the nucleus, lying behind the stars 
to the E. This geometry is consistent with a stellar bar with leading dust lanes. Greve et 
al. (2002), based on extinction and polarization arguments, suggested that M82 is observed 
from below. 



5.3.1. The presence of a stellar bar. 

The total H92a line distribution extends over ~ 35" along the major axis of M82 (Fig 5), 
covering the inner plateau and also the inner edges of the outer plateau. Since the resulting 
velocity field provides the average velocity value of the gas at each position, we have used 
the data cube to search for the steepest velocity gradient in the center of M82. The steepest 
velocity gradient is determined from terminal velocities (~ 36 km s^^ arcsec^^) and is found 
along a PA of ~ 57°. Only two regions, separated by ~ 7'.'8, define this velocity gradient: one 
of these regions is Ca and the other is a region located 0'.'4 from Wlc. The center defined by 
these two regions is located 1''7 E from the 2.2 //m peak determined by Dietz et al. (1986). 
The velocity gradient could be alternatively defined between the region Ca and a region 
located near the center of Ela; in this case a velocity gradient of ~ 27 km s~^ arcsec"^, 
at a PA of ~ 68° is obtained. This PA agrees with the major axis of M82 determined on 
large scale. The position at mid-distance between these two regions (separated by ~ 10") 
has an offset of = — 1"5, As = 0"4 relative to the 2.2 //m peak (Dietz et al. 1986). This 
offset corresponds approximately to the distance between the 2.2 /xm peak and the center 
of the inner plateau observed in the 2.2 jim. profile (see Fig. 2 of Larkin et al. 1994). This 
asymmetry, with a truncation of the inner bar on the E side of the 2.2 /im profile, has already 
been noticed by Larkin et al. (1994) and could reveal the presence of ionized gas in two lobes; 
this suggestion is supported by our H92q; observations since the regions Ca and Ela appear 
as the only significant features in the channels at 313 and 284 km s~^. Although different 
velocity gradients are obtained using the regions near Ela and Wlc, both regions are located 
on the ridge of extinction determined by Larkin et al. (1994). In any case, the center position 
defined by these two regions is very close to the location of maximum extinction along this 
ridge. Dust lanes are expected to be present at the intersection of xl and x2 orbits. This 
result suggests the presence of elongated x2 orbits. Assuming that the PA of the xl orbits 
corresponds to the outer plateau is ~ 70° (Telesco et al. 1991), the major axis of the x2 
orbits would be close to an end-on view. 

The position-velocity {p-v) diagram, using the integrated fiux density along the minor 
axis direction by taking shces at a PA of 70°, is shown in Figure 10 (top). This PA was used 
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in order to compare with results from [Ne II] and Br7 observations (Achtermann and Lacy 
1995; Larkin et al. 1994). In the H92q; hne image there is emission at "forbidden" velocities 
(near —I'.'S at 215 km s~^), inconsistent with purely circular motions. If non-circular motions 
are present, additional velocity components are observed. This feature is also observed in 
the Br'-f line images (Fig. 5 of Larkin et al. 1994). In [Ne II], there is also emission at 
"forbidden" velocities but this is distributed in a symmetrical fashion about 200 km 
around the 2.2 fxm peak. The agreement of the H92q; p-v image is closer to the Br^y (Larkin 
et al. 1994) p-v images than to the [Ne II] p-v image. 

Figure 10 (bottom) shows the p-v image constructed along the major axis using the PA 
of 65°. In this case the profiles are in general relatively narrower when they are compared 
with those in Fig. 10 (top); the binning along the minor axis is over a smaller region of ~ 1''5 
as compared with 9". However, the spectra remain relatively broad except near 8" on the W 
side and 10" on the E side. These two regions are located at the transition zones between 
the inner and outer 2.2 fim plateaus. The asymmetry observed at 2" angular resolution, 
relative to the position of the 2.2 //m peak, is still present for this PA=65°, suggesting that 
the asymmetry is very likely to arise in the bar itself. In this case, there is no emission 
at "forbidden" velocities, indicating that the orbits must be confined to a relatively thin 
plane not observed edge-on. By comparing with optical measurements (McKeith et al. 1993; 
Greve et al. 2002), we observe systematic effects: in the inner part (< 6") the terminal H92q; 
velocities are close to the optical velocities only on the W side. Neither the gas nor the stars, 
as observed in the optical, have velocities near ~ 325 km s~^ from to -\-6". On the other 
hand, there are no optical velocities near ~ 200 km from —1" to —6". Beyond ±11", on 
the outer 2.2 /xm plateau, the optical velocities are systematically closer to the low velocity 
component of the II92a profiles. These results would indicate that for the regions located 
in the outer plateau, the gas as well as the stellar emission observed in the optical comes 
preferentially from the leading edge of the bar in the xl orbits (see Fig. 2 of Greve et al. 
2002). 

The symmetry observed in the terminal velocities of the H92« line would indicate that 
the component of the ionized gas that we observe is not only in the xl orbits but also in x2 
orbits (< ±3") and/or with some "spraying" (at radial distances from 3" to 6"). The result 
of gas shocking in the intersection of xl and x2 orbits is named by some authors as "spray" 
orbits (e.g. Downes et al. 1996). The encounter of gas moving in these two orbits breaks 
the flow along the xl orbit into a spray that traverses interior of the bar until it reaches the 
far side of the bar. A schematic diagram for these spray orbits is shown by Downes et al. 
(1996). In this context it is difficult to explain why the optical velocities are close to the 
H92q; terminal velocities only on the W side. One possibility would be that the gas observed 
in the optical does not spray on the E side and the gas emitting the H92q; line has orbits 
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with radii which extend relatively closer to the nucleus. An alternate possibility would be 
that the x2 orbits are in a plane closer to 81° inclination than the xl orbits, the far side 
of the x2 orbits emitting in the optical being more severely affected by the extinction from 
dust in the plane of the xl orbits. 

There is a thin transition edge between the inner and outer plateaus not only observed 
in the p-v images made along the major axis (Fig. 10) but also in the velocity dispersion 
image (Fig. 8b). In Figure 8b an ellipse (fitted through the regions with local minima) is 
shown. This ellipse is centered at the same position as the ring found by Achtermann and 
Lacy (1995), with a major axis of 17" at a PA of 72° and the same axial ratio of ~3.5. 
This result agrees with the 2.2 /^m light distribution for the inner plateau, consistent with a 
picture in which the inner and outer plateaus correspond to different families of orbits. 

5.3.2. The inner ring. 

A description in terms of elongated xl and x2 orbits is not the only model that can 
explain the observed kinematics in M82. The observed velocities distribution is not consistent 
with a rotating disk unless it consists of a set of rings with different tilt angles which may 
be contracting. Different methods to fit a ring in the center of M82 were used; the results 
are fisted in Table 6. Row 1 of Table 6 lists the results of Achtermann and Lacy (1995), 
column 1 lists the central position of the ring, column 2 the position angle of the major axis, 
columns 3 and 4 the angular size of the major and minor axis, column 5 the systemic velocity, 
column 6 the rotation velocity, column 7 the expansion velocity, column 8 the dispersion 
about the model and column 9 the spectral line used to obtain the parameters of the ring. 
Achtermann and Lacy (1995) fixed the geometry of the ring based on the physical appearance 
of the distribution of the [Ne II] line emission and then fitted the velocity pattern along that 
ring, assuming that the ring lies in the plane of the galaxy, centered on the kinematical 
center. Row 2 is the result obtained from the H92q; velocity pattern (Fig. 8a), using the 
same morphological parameters listed in row 1. Following the same procedure for the H92q; 
line emission, the ellipse fitted by Achtermann and Lacy (1995) would need to be shifted 
~ O'.'S to the N (Fig. 4) to coincide with the ridge of emission. In row 3, this shift of 0''3 has 
been applied. In row 4, the position of the center is a free parameter. Finally, row 5 lists 
the parameters of the ring (shown in Fig. 4) obtained without imposing any morphological 
constraints; this procedure gives the solution with the lowest residual dispersion. 

Based on the spatial distribution of the line emission, the existence of the ring is not 
obvious in any of the models presented in Table 6. According to the ring model of Achter- 
mann and Lacy (1995) and the one described by row 5 (Table 6), the rotation velocity Vrot 
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has a first maximum at a radius of ~ 5" — 6". In the model of Achtermann and Lacy (1995) 
the ring surrounds a region devoid of ionized gas, in contrast with the model presented here 
(the best fit to the velocity pattern). Both models require not only a contraction velocity 
but also that the inclination angle of the rings is different than the inclination of the plane 
of the outer part of M82 (i=81°). The major axis of these rings is very close to the PA 
determined from the large scale morphology (68°). 

In the velocity dispersion image (Fig. 8b), a well defined band at a PA of ~ 150° is 
observed with systematically narrower line profiles. The length of this band is not restricted 
to the extent of the minor axis of the ring with axial ratio of ~ 3.5 (Achtermann and Lacy 
1995); in addition this band covers the extent of the velocity pattern along that direction. 
On both sides of this band, there is a ridge in which the broadest profiles are observed. 
These regions of large velocity dispersions are explained by the terminal velocities with 
|Vt — Vsys\ — 120 km s^^, along with a second velocity component in the spectra with 
velocities closer to Vsys, |Vt — Vsys\ — 50 km s~^. Figure 11 shows the p-v diagram along 
the E ridge. The regions where Vt > +300 km s~^ are not restricted to the range < ±1''7, 
which corresponds to the regions associated with the ring of axial ratio of ~3.5 (see Table 
6). Outside the range —5" to +3", i.e. beyond the ring of axial ratio 2 (see Table 6), the 
bright component remains confined near ~ 250 km s~^. We conclude that, if gas in x2 orbits 
and/or a fast rotating nuclear ring (spray between the xl and x2 orbits) is responsible for 
Vt > +300 km s~^, it cannot be confined in a thin plane which would give this axial ratio of 
~ 3.5. 



5.3.3. Outflow in the halo. 

Fig. 11 shows that for offsets larger than —5" there is emission at two different velocities, 
one component emerging near ~ 150 km s~^ and the second one near ~ 350 km s~^. At 
~ —10" the emission is detected in the range ~ 250 km and ~ 30 — 100 km s~^. These 
results agree with the optical measurements (McKeith et al. 1993) in the lines of [SII], [Nil], 
Ha. A more detailed view of the radial velocities of the outflow, observed in Ha, are shown 
in Figure 10 of Shopbcll et al. (1998). The H92« observations reveal the base of the outflow 
into the halo. On the N side of M82 there are two velocity components: the ionized gas on 
the far side of the cone (component I) has radial velocities ~ 250 km s^^ and on the near 
side of the cone (component II) the gas has radial velocities in the range ~ 30 to 100 km s~^. 
On the S side (~ +5"), the velocities are confined to ~ 200 km s^^ (component IV), again 
in agreement with the optical measurements (Shopbell et al. 1998). The geometry of the 
outflow and description of these four components are shown in Figure 5 of McKeith et al. 
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(1995). 

Due to the high inchnation of M82, it is difficult to distinguish between an ad-hoc 
model of circular motions with radial contraction vs elliptical xl/x2 orbits. The ring models 
indicate contraction instead of expansion velocities (Table 6). The presence of different 
velocity components in the H92q; fine emission along a PA of 150° indicates that the ionized 
gas near the nucleus has already been perturbed by the outflow. A barred potential seems to 
provide the mechanism to bring gas from the outer parts of the disk of M82 into the central 
part. In this way, the gas flowing into the center of M82 would be consumed in the starburst 
with some gas outflow into the halo. 

5.4. Existence of an AGN 

The presence of an AGN in the center of M82 has been proposed by Matsumoto & 
Tsuru (1999) based on X-ray data. These authors suggest that there is an under- luminous 
AGN in the center of M82, comparable to Sgr A*. In the continuum at 8.3 GHz (Fig. 2a), 
wc observe a source located near the 2.2 fim peak (Lester ct al. 1990). In the continuum 
at 43 GHz, wc do not detect this source above the 3arms level. Muxlow et al. (1994), using 
MERLIN observations at 5 GHz with angular resolution of 50 mas, report a compact source 
(43.55-1-60.0) that was identified as a young supernova remnant with shell-hke structure. 
Thus, no clear evidence has been found for the existence of a compact radio source that 
can be associated with an AGN in the center of M82 within our 8.3 GHz Sarms limit of 
0.12 mJy beam~^. 

6. CONCLUSIONS. 

We have presented high angular resolution (< 1") VLA observations of radio continuum 
emission from M82 at 8.3 and 43 GHz and of the H92q; and H53q; RRLs. In the radio 
continuum images, we identified 19 new compact sources at 8.3 GHz and five at 43 GHz. 
Of these newly identified sources, four are SNRs, four are HII regions and the nature of the 
remaining 13 is uncertain. Three SNRs show the presence of H92a; RRL emission which can 
be interpreted as line emission produced by HII regions along the line of sight. No compact 
source is observed near the position of the derived kinematic center. Thus, in agreement 
with Muxlow et al. (1994), we find that there is no compeUing evidence for the existence of 
an AGN at the center of M82. 

We have modeled the line and continuum emission using a collection of HII regions. We 
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have considered models with a single density as well as models with two density components. 
The H53q; RRL mainly traces the high-density ionized gas. Models with multiple HII regions 
were used to reproduce the observations in the continuum at 5, 8.3 and 43 GHz as well as the 
H92a and H53a RRL emission from four separate Hll complexes. The two-density model is 
considered to be more appropriate to model the line and continuum emission since compact 
continuum sources have been identified already as Hll regions. However, as the two-density 
model is based on previous identification of thermal compact sources, the emission from 
region El (no compact thermal sources have been identified) has been modeled using only 
the single-density model. In general the low-density component is characterized by HE 
regions with sizes of ~0.8 pc and an average electron density of ~ 5 x 10^ cm~^ and the 
high-density component is characterized by HII regions with typical sizes of ~ 0.15 pc and 
an average electron density of ~ 3.5 x 10'* cm^'^. From the two-density model, we can infer 
that the starburst in the region Wl is the youngest. The derived mass of ionized gas in the 
central region of M82 is ~ 2 x 10^ M©, 30 times lower than for Arp 220 (Anantharamaiah 
et al. 2000). Based on the models, it has been inferred that M82 has a star formation rate 
of ~ 3 Mq year~^, a factor ~ 100 lower compared with the merging system Arp 220. For 
a more accurate determination of the properties of the ionized gas over a complete density 
range, high-resolution observations of RRLs at additional higher and lower frequencies are 
required. 

The H92q; line emission extends over ~ 35". The steep velocity gradient, measured at 
PA of 68°, is 26 km s~^ arcsec"^ in the inner ~100 pc of M82. At a distance of ~80 pc 
from the kinematic center, the gas has a maximum radial velocity of ~150 km s~^, implying 
an enclosed total dynamical mass of ~10^ M0. From the p-v diagram of the H92a line, 
we observe deviations from circular motions. The observed velocity pattern cannot be due 
to pure circular motions; apparently the orbits in the inner parts of M82 arc more face-on 
than in the outer parts. In the velocity dispersion image, the H92a line traces the base of 
a large scale outflow into the halo. This outflow is observed only on the N side. We find 
the PA of the axis of the outflow to be 150°. Along this axis, the widths of the profiles 
are narrow, suggesting coUimated flow with small opening angle in the regions within the 
disk. Within the disk, the steepest gradient (~ 36 km s~^ arcsec"^), is measured from the 
terminal velocities. This gradient is observed at a PA of ~ 57°, in a direction perpendicular 
to the axis of the outflow into the halo. This direction is not parallel to the PA of the large 
scale major axis of M82 (~ 65°) nor the PA of the 2.2 /im bar (~ 70°). 

The best fit with an ad hoc model of circular orbits leads to a contracting ring that 
has its kinematic center shifted ~ 1" W of the 2.2 //m peak and located approximately in 
the middle of the 2.2 /xm inner plateau. The nodes characterizing this ring coincide with 
the dust lanes; this coincidence could represent additional evidence of gas in x2 orbits of a 
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bar potential. Due to the presence of the outflow, the kinematics observed in the H92ci; hne 
provide possible evidence for the presence of x2 orbits originating from a bar potential. The 
inner and outer 2.2 ^m plateaus correspond to two different families of orbits; the transition 
zone between the two plateaus has a signature in the kinematics in the form of a reverse in 
the velocity gradient that has been observed near the major axis. 

The National Radio Astronomy Observatory is a facility of the National Science Foun- 
dation operated under cooperative agreement by Associated Universities, Inc. We thank J. 
M. Torrelles for making available his H92« data. We also thank Yolanda Gomez and anony- 
mous referee for their very useful comments. CR acknowledges the support from UNAM 
and CONACyT, Mexico. 
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Tabic 1. PARAMETERS OF THE VLA OBSERVATIONS OF M82 



Parameter 


H92a Line 


H53q; Line 


Right ascension (B1950) 


09'^51'"42?38 


09'*51'"42?51 


Declination (B1950) 


69°54'59^'2 


69°55'00^'0 


Total observing duration (hr) 


32 


13 


Bandwidth (MHz) 


25 


50x2 


Number of spectral channels 


31 


15x2 


Center YrcI (km s~^) 


200 


160 


Velocity coverage (km s"-*^) 


700 


350x2 


Velocity resolution (km s~^) 


56 


44 


Amplitude calibrator 


3C286 


3C286 


Phase calibrator 


1044+719 


0954+658 


Bandpass calibrator 


3C48 


1226+023 


RMS line noise per channel (mJy/beam^ ) 


0.07 


1.3-2.0 


RMS, continuum (mJy/beam'^ ) 


0.04 


0.5 



^Synthesized beam size of 2". 
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Table 2. CONTINUUM FLUX DENSITIES, COMPACT FEATURES. 



8.3 GHz 43 GHz 



Source ID 


Coordinates 


c- 2 

Size 


S8.3 


Coordinates 


c- 2 
Size 


S43 


Spectral index 


RRL 


Type 






( ) 


(mjy) 




{ ) 


/ T \ 

(mJy) 


O.0-4.5 LrXlZ 


Low freq. 


feature 




(1) 


(2) 


(3) 


(4) 


fr:\ 

(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


37.54+53.2 


37.54+53.2" 


0.59 


0.43+0.03 






<1.5 


<0.74 








38.76+53.4 


38.75+53. 5*" 


0.5s 


0.24+0.07 






<1.5 


<1.1 


0.54±0.27<^ 




HH.'^ 


39.10+57.3 


39.11+57.3*' 




3.81+0.04 






<1.5 


< —0.55 


—0.53+0.09'' 




SNR-'' 


39.29+54.2 


39.27+54.1*' 


0.4 


3.6+0.1 


39.25+54.1*' 


0.5 


3.5+0.3 


-0.02+0.02 


1.64+0.16'' 




HII-'^ 


39.40+56.1 


39.43+55.9' 


0.5 


1.9+0.1 






<1.5 


< -0.14 


-1.04+0.22'* 




SNR-'' 


S9 64-1-5'^ 4 


39 65-^53 1^ 




1 0+0 05 






<1.5 


<0.24 


-0.20+0.05'= 




SNR-'' 


39.68+55.6 


39.67-1-55.3'' 


0.4 


3.9+0.2 


39.66+55.4*' 


0.5 


3.6+0.5 


-0.05+0.09 


1.03+0.09'* 




HII.'^ 


SQ 77-1-56 Q 


39 78-1-56 9*' 


0.4 


7+0 05 






<1.5 


<0.45 


—0.50+0.06'= 




SNR-'' 


39.92+55.9 


39.92+55.9" 


0.7 


2.7+0.1 






<1.5 


< -0.35 






SNR?== 


40.10+55.0 


40.15+54.5'' 


0.7 


0.6+0.05 






<1.5 


<0.54 








40.32+55.1 


40.33+55.1'' 


0.4 


0.8+0.04 






<1.5 


<0.37 


-0.23+0.21=* 




SNR-'' 


40.49+57.4 


40.49+57.4" 


1.0 


2.7+0.1 






<1.5 


< -0.35 






SNR?'= 


40.62+56.0 


40.63+56.0*' 


0.5 


1.7+0.1 






<1.5 


< -0.01 


-0.72+0.25=* 




SNR-' 


40.66+55.2 


40.67+55.1*' 


0.4 


5.6+0.1 


40.65+55.1*' 


0.5 


2.7+0.1 


-0.43+0.05 


-0.54+0.08=* 




SNR-f 


40.95-1-58.8 


40.93+58.7*' 


0.5 


4.5+0.4 


40.92+58.7*' 


0.5 


3.9+0.2 


-0.09+0.06 


0.44+0.13=* 


W2h 


HH-' 


40.96+57.9 


40.96+57.8*' 




2.3+0.2 


40.96+57.6*' 


0.5 


1.8+0.1 


-0.15+0.05 


>1.2+ 






41.17+56.2 


41.17+56.1*' 


0.4 


5.4+0.2 


41.16+56.2*' 


0.7 


4.9+0.4 


-0.06+0.03 


0.87+0.14'' 


W2f 


HII-' 


41.29+59.7 


41.33+59.2*' 


0.7 


3.0+0.2 






<1.5 


< -0.41 


-0.47+0.09=* 


W2e 


SNR-f 


41.62+59.9 


41.62+59.9" 


0.5 


1.6+0.1 






<1.5 


< —0.04 








41.95+57.5 


41.95+57.4*' 




13.6+0.4 


41.95+57.5*' 


0.7 


5.0+0.1 


-0.60+0.02 


-0.80+0.05= 




SNR?-'' 


42 08-1-58 4 


42 11+58 3*' 


0.4 


4.0+0.2 


42 09+58 3*' 


0.5 


3.6+0.2 


-0.06+0.04 


1.32+0.17=* 




HH-' 


42.21+59.0 


42.20+59.0*' 


0.5 


6.7+0.1 


42.21+59.0*' 


0.7 


6.5+0.3 


—0.02+0.03 


1.16+0.13=* 


Wlf 


HH-' 


42.67+55.6 


42.67+55.5*' 




0.92+0.04 






<1.5 


<0.29 


-1.3+0.2=* 




SNR-f 


42 69+58 2 f 
















] 1.04+0.13 


Wlc 


HIl/ 


42.56+58.0 1 


42.64+57.9*' 


~ 1.4 


32.4+0.7 


42.64+57.9*' 


~ 1.2 


16.2+0.8 


-0.41+0.2 


1 0.88+0.14 


Wlc 


HII-' 


42.48+58.4 [ 
















J >1.2 


Wlc 


HH-' 


43.00+59.0 


43.02+58.9*' 


0.4 


3.9+0.2 






<1.5 


< -0.56 




Wla 


SNR?= 


43.18+58.3 


43.17+58.3*' 


u 


5.5+0.1 






<1.5 


< -0.77 


-0.44+0.08=* 




SNR-f 


43.21+61.3 


43.21+61.3" 


0.8 


4.3+0.1 


43.25+61.4" 


1.3 


4.0+0.3 


-0.04+0.05 






Hn?= 


43.31+59.2 


43.29+59.0*' 


u 


8.0+0.2 


43.30+59.2*' 


0.4 


2.0+0.2 


-0.83+0.06 


-0.65+0.07=* 




SNR-f 


43.39+62.6 


43.39+62.6" 


0.7 


7.0+0.2 


43.39+62.7" 


0.8 


3.9+0.2 


-0.34+0.14 




Cd 




43.50+61.2 


43.50+61.2" 


0.5 


2.2+0.1 






<1.5 


< -0.23 








43.55+60.0 


43.57+59.8*' 


0.4 


0.77+0.05 






<1.5 


<0.40 








43.65+57.7 


43.65+57.7" 


0.7 


2.3+0.1 






<1.5 


< -0.25 








44.01+59.6 


44.02+59.5*' 


0.4 


19.5+0.2 


44.01+59.5*' 


0.4 


4.5+0.4 


-0.87+0.05 


-0.38+0.06=* 




SNR-f 


44.08+63.1 


44.08+63.1" 


0.8 


4.1+0.1 






<1.5 


< -0.60 








44.11+64.3 


44.11+64.3" 


0.5 


2.1+0.1 






<1.5 


< -0.20 






SNR?'= 


44.17+64.4 






< 0.24 


44.17+64.4" 


0.8 


4.1+0.2 


> 1.74 






HIP 


44.29+59.3 


44.28+59.3*' 


u 


2.3+0.1 






<1.5 


< -0.25 


-0.72+0.12=* 




SNR-f 


44.36+57.8 


44.36+57.8" 


0.4 


1.5+0.1 






<1.5 


<0.0 








44.43+62.5 






< 0.24 


44.43+62.5" 


0.9 


5.5+0.2 


> 1.91 




Elb 


HIP 
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Table 2 — Continued 



8.3 GHz 43 GHz 



Source ID^ 


Coordinates^ 


Size2 


Ss.s 


Coordinates-'^ 


Size2 


S43 


Spectral index 


RRL 


Type 






(") 


(mjy) 




{") 


(mJy) 


8.3-43 GHz 


Low freq. 


feature 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


44.0U-r00.o 


A A c^n_L^;c; oa 
4l:4.0U-rt)0.o 


u.o 


1 /I-Un 1 

1.4l:ltU.l 






<i.O 


<0.0 








44.52+58.1 


44.52+58.0*' 


0.4 


2.1+0.05 






<1.5 


< -0.20 


-0.15±0.17<^ 




SNR-f 


44.84+61.8 


44.84+61.8" 


u 


0.95+0.05 






<1.5 


<0.27 








44.91+61.1 


44.89+61.2* 


u 


1.2+0.1 






<1.5 


<0.13 


-0.45+0.18"* 




SNR^ 


44.93+63.9 


44.92+63.7*' 


0.7 


3.6+0.1 






<1.5 


< -0.52 


>1.1'* 


E2f 


HH-^ 


45.17+61.2 


45.19+61.1* 


0.4 


6.7+0.1 






<1.5 


< -0.89 


-0.52+0.07'' 




SNR/ 


45.38+60.3 


45.38+60.3'^ 


0.4 


1.1+0.05 






<1.5 


<0.19 








45.44+67.3 


45.44+67.3' 


0.4 


1.47+0.06 






<1.5 


<0.0 








45.63+66.9 


15.63+66.9° 


0.8 


5.8+0.2 


45.61+67.1° 


0.9 


5.5+0.4 


-0.04 




E2c 


Hn= 


45.70+62.9 


45.66+62.8* 


0.4 


0.78+0.05 






<1.5 


<0.39 








45.79+65.2 


45.74+65.5* 


0.5 


5.9+0.2 






<1.5 


< -0.81 


-0.55+0.13'' 


E2b 


SNRf 


45.91+63.8 


45.88+63.7* 


0.4 


1.28+0.05 






<1.5 


<0.01 


-0.38+0.11'* 




SNR/ 


45.93+74.3 


45.93+74.3° 


1.0 


1.4+0.1 






<1.5 


<0.01 








46.17+67.6 


46.17+67.6'' 


0.8 


4.6+0.1 


46.21+67.8* 


0.8 


4.5+0.3 


-0.02+0.04 


0.66+0.16'' 


E2a 


HH-'^ 


46.33+66.2 


46.33+66.2° 


u 


0.93+0.05 






<1.5 


<0.28 








46.52+63.8 


46.53+63.9* 


0.4 


2.1+0.1 






<1.5 


< -0.20 


-0.35+0.11"* 




SNRf 


46.56+73.8 


46.57+73.7* 


u 


0.76+0.05 






<1.5 


<0.40 


-0.60+0.05'= 




SNRJ' 


46.74+69.7 


46.74+69.7" 


1.4 


4.3+0.2 






<1.5 


< -0.63 






SNR?'= 


46.75+67.0 


46.69+66.9* 


0.4 


1.8+0.1 






<1.5 


<0.11 


-0.80+0.05= 




SNR-*^ 


47.11+66.3 


47.11+66.3° 


0.5 


1.12+0.05 






<1.5 


<0.17 








47.37+68.0 


47.37+67.8* 


0.4 


0.7+0.05 






<1.5 


<0.45 


-0.60+0.07= 




SNR/ 



^Relative to (B1950) 09''51'"+ 69° 54'. 

^Deconvolved angular size from the 0.6" HPBW. With this resolution the sources which appear unresolved are indicated by "u". 
''Sources newly identified in this paper. 

''Prom Kronbcrg ct al. (1985); Bartel et al. (1987); Huang ct al. (1994); Muxlow et al. (1994); Wills et al. (1997); McDonald et al. (2002). 
"^Spectral index (5 oc i/") between 1.4, 5 and 8.4 GHz (Wills et al. 1997). 
''Spectral index (5 (X i/") between 5 and 15 GHz (McDonald et al. 2002). 
^Source type obtained in this paper based on the spectral index in column 8. 
'Source type given by McDonald et al. (2002). 



Table 3. LINE PARAMETERS AT LOW ANGULAR RESOLUTION USED IN THE MODELS. 



H92a H53a 



Region 


Peak Flux 
(mJy) 


^HEL 
(km/s) 


AV 
(km/s) 


Sc8 
(mJy) 


Peak Flux 
(mJy) 


^HEL 
(km/s) 


AV 
(km/s) 


So43 
(mJy) 


(10-3) 


Sc6 
(mJy) 


E2 


2.2±0.1 


315±2 


106±5 


106±10 


7.5±1.5 


305±5 


55±12 


46±17 


1.3 


152 


El 


l.OiO.l 


269±3 


130±8 


72±10 


6.7±1.5 


264d=5 


50±13 


28±15 




93 


Wl 


6.4±0.1 


119±1 


116±3 


292±15 


9.8±2.0 


113±12 


120±29 


157±20 


3.2 


403 


W2 


4.2±0.1 


100±2 


112±4 


169±15 


10.2±2.0 


66±9 


100±22 


101±19 


6.8 


230 



*Note that the integrated line flux density for eaeh region is lower than that given in Table 4 (column 2) since the regions 
for which we have integrated here are defined by the lowest contour in the H53a integrated line emission image. 

''Area filling factor computed using Table 4 from McDonald et al. (2002). 
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Table 4. PARAMETERS FOR REGIONS OBSERVED IN H92a RRL 





Low angular resolution 






Intermediate resolution 






Region 


/SidV 




Features 


Coordinates 


/SidV 


Size= 


Peak Flux 


V HEL 


/\\/d 

FWHM 




(mjy km s~^) 


(xlO-2) 






(mJy km s 


-1) („) 


(mJy) 


(^Km/s ) 






(2) 


(3) 


\^) 




(6) 


(7) 


(8) 


V') 




E2 


604 


1.7±0.1 


E2a 


46.20+67.7 


81 


1.3 


0.57+0.03 


316+16 


121+37 








E2b 


45.79+65.3 


65 


1.4 


0.82+0.05 


315+5 


48+11 








E2c 


45.54+66.3 


165 


1.6 


1.54+0.14 


320+8 


84+19 








E2d 


45.43+65.6 


96 


0.8 


0.82+0.05 


311+5 


95+11 








E2e 


45.32+64.2 


69 


0.9 


0.61+0.05 


314+5 


90+11 








E2f 


45.05+63.5 


32 


0.4 


0.26+0.06 


302+7 


102+16 








E2a+E2b+...+E2f 




508 










El 


300 


1.2±0.1 


Ela 


44.39+61.8 


98 


1.4 


0.87+0.07 


264+8 


90+18 








Elb 


44.10+62.7 


56 


1.6 


0.56+0.05 


321+10 


76+24 








Elc 


44.02+60.3 


39 


1.4 


0.59+0.05 


231+10 


24+24 








Ela+Elb+Elc 




193 










C 


320 


l.OitO.l 


Ca 


43.84+62.4 


70 


0.9 


0.63+0.05 


300+5 


88+11 








Cb 


43.77+63.2 


137 


1.2 


0.61+0.05 


286+5 


205+11 








Cc 


43.79+63.9 


16 


0.6 


0.30+0.05 


470+5 


52+11 








C(l 


13.39+62.3 


75 


1.3 


1.00+0.05 


239+4 


44+9 








Ca-Cb+Cc-Cd 




2l)iS 










Wl 


1067 


1.6±0.1 


Wla 


42.95+59.2 


85 


1.9 


0.72+0.10 


165+8 


96+4 








Wlb 


42.62+60.3 


196 


2.1 


1.85+0.10 


148+3 


83+8 








Wlc 


42.61+58.0 


291 


1.7 


2.88+0.07 


117+2 


77+3 








Wld 


42.32+61.5 


57 


1.1 


0.56+0.10 


191+8 


78+4 








Wle 


42.21+59.9 


147 


1.3 


1.05+0.07 


125+2 


120+3 








Wlf 


42.19+58.8 


177 


1.6 


2.01+.10 


100+3 


61+8 








Wla+Wlb+...+Wlf 




953 










W2 


878 


2.0±0.1 


W2a 


41.75+59.0 


113 


1.4 


0.76+0.07 


114+12 


128+32 








W2b 


41.62+57.8 


141 


u 


1.43+0.07 


79+12 


69+32 








W2c 


41.51+58.6 


81 


1.2 


0.68+0.07 


113+3 


96+6 








W2d 


41.47+56.6 


69 


1.6 


0.80+0.08 


76+5 


57+12 








W2e 


41.27+59.2 


65 


1.1 


0.69+0.11 


135+5 


68+10 








W2f 


41.12+56.4 


113 


1.8 


1.19+0.08 


86+5 


69+12 








W2g 


41.07+61.0 


36 


1.3 


0.23+0.07 


84+3 


133+6 








W21i 


40.97+58.5 


112 


1.4 


1.14+0.09 


129+4 


73+9 








W2a+W2b+...+W2h 




730 











^The So and Sj^ values are listed in Table 3. 
''Relative to (B1950) 09''51™+ 69° 54'. 

'^Deconvolved angular size from 0"9 synthesized beam, u is given when the source is unresolved. 
■^Deconvolved line width from 56 km ~^ , the spectral resolution achieved in the H92a line. 
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Table 5. MODELS BASED ON THE H92ci; AND H53q; LINES 



Parameter 


Region E2 


Region El 


Region Wl 


Region W2 






SDd 







Te X 10^ (K) 7±2 

ne X 103(cm-3) .... 9±8 

Size (pc) 0.81±0.59 

Nhii X 10^ 2.5±1.0 

P 0.29±0.28 

Log Nlj; 52.36±0.19 

SFRb (M0 year-l) 0.37±0.13 

rc (8.3 GHz) 0.27±0.26 

tl (8.3 GHz) -0.13 ±0.10 

b„ (H92a) 0.97±0.02 

/?„ (H92a) -20±15 

a" 0.95±0.20 

Mhii X 10-* (Mq) . . 2.5±2.6 



7.5±2.5 
6±1 



40 ±20 
5±4 



0.14±0.06 
0.86±64 



0.56±0.36 
2.4±1.1 



0.0013 
0.29±0.27 



51.64±0.32 
52.33±0.21 



().1()±().07 
0.45±0.19 



0.032±0.023 
0.14±0.13 



-0.050±0.048 
-0.095±0.068 



0.92±0.03 
0.96±0.02 



-23±8.1 
-22±13 



1.0 ±0.24 



0.041±0.008 
3.0±2.6 



7.5±2.5 
26±25 
0.75±0.65 

1.7±0.88 
0.27±0.26 
52.16±0.25 
0.21±0.11 
0.51±0.50 
-0.16 ±0.13 
0.97±0.03 

-20±16 
1.00±0.30 

1.5±1.4 



6±0.5 
4±3 
0.72±0.34 

5.8±1.1 
0.31±0.24 
52.75±0.09 
1.64±0.13 
0.10±0.08 
-0.09 ± 0.05 
0.97±0.02 

-19±9 
0.50±0.03 
6.6±4.5 



r 7.5±2.5 
\ 5±0.5 

f 30±10 
\ 4.0±2.5 

■ 0.15±0.04 
0.73±33 

f 1.6±0.7 
\ 4.6±0.3 

0.0032 
0.30±0.23 

52.15±0.2 
52.67±0.()3 

/ 0.29±().13 
\ 0.86±0.05 

f 0.015±0.007 
\ 0.09±0.07 

-0.032±0.027 
-O.()85±0.042 



{ 



0.91±0.025 
0.97±0.014 



-26±5 
-20±8 



0.52 ± 0.026 

/ 0.16±0.02 
\ 6.1±4.0 



7.5±2.5 
15±13 
0.52±0.36 

6.8±3.8 
0.37±0.35 
52.75±0.27 
0.63±0.51 
0.35±0.32 
-0.14 ±0.11 
0.98±0.02 
-18±12 
0.91±0.5 
5.6±5.1 



7.5±2.5 
6±1 



30±10 
5.5±4.0 



0.15±0.04 
0.69±0.38 



2.4±1.0 
4.6±2.6 



0.0068 
0.40±0.37 



52.34±0.2 

52.58±0.27 



0.45±().19 
0.86±0.48 



().015±0.007 
0.14±().12 



-().()51±0.047 
-().()95±0.068 



{ 



0.91±0.025 
0.97±0.02 



-26±5 
-20±11 



{ 

0.62 ±0.19 

I 0.25±0.03 
\ 6.1±6.4 



*Area covering factor. 

''Star formation rate using the formulae given by Anantharamaiah et al. (2000). 
'^Spectral index between 8.3 and 43 GHz. 

''Results from single-density models are listed in columns SD and results from two-density models are listed in columns labeled TD. 
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Table 6. MODELS FOR A NUCLEAR RING. 



Ring center ' 


PA 


Maj 


Min >> 


Vsys 


Vrot 


Vexp 


"'fit 


Line 




(deg.) 


(arcsec) 




(km E 








9''51"'43f40, 69°55'0;'l'' 


70 


11.4 


3.3 


200 


112 


-16 




[Ne 11]'^ 


9''51'"43f40, 69°55'0;'l'' 


70 


11.4 


3.3 


193 


94 


-13.5 


7.47 


H92Q 


9''51™43?40, 69°55'0'.'4'' 


70 


11.4 


3.3 


196 


96 


-15 


7.27 


H92o 


9'"51'"43fl9, 69°55'0!'6 


70 


11.4 


3.3 


191 


102 


-13 


7.35 


H92a 


9''51'"43f28, 69°55'0f'6 


66 


11.4 


5.6 


195 


103 


-14 


4.82 


H92a 



^Results from Achtermann and Lacy (1995). 

''The fixed parameters are given in bold face and free parameters estimated from the fit in normal face. 
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Fig. 1. — a) Radio continuum images at 8.3 GHz made with the VLA in the B and C arrays. 
Contour levels are drawn at -3, 3, 269 times the rms noise of 0.07 mJy beam~^ at steps 
of 2^/^. The gray scale covers the range from to 70 mJy beam"^. b) Radio continuum 
images at 43 GHz made in the C array of the VLA. Contour levels are -3, 3, 4, 6, 8, 10 times 
0.4 mJy beam~^, the rms noise in the central regions. The gray scale image covers the range 
from 3 mJy beam~^ to 53 mJy bcam~^. In both images, the radio continuum made with 
high angular resolution of 0''6 (contours) is superposed on the radio continuum image with 
low angular rcsohition of 2" (gray scale). The square shows the position of the 2.2 /xm peak 
(Lester et al. 1990) and the circle marks the position of the SNR G41. 95+57.5. 




Fig. 2. — Radio continuum images after subtracting the background local emission, (a) 
8.3 GHz image (B and C arrays). Contour levels are -3, 3, 6, 12, 24, 48, 96 times the rms 
noise of 0.04 mJy beam~^. (b) 43 GHz image (C array). Contour levels are drawn at -3, 3, 
4, 6, 8, 10 times the rms noise of 0.4 mJy beam~^. Both images have angular resolution of 
0''6. The crosses (+) indicate the position of SNRs, while the rotated crosses (X) indicate 
the HH regions (McDonald et al. 2002). The empty squares show the position of sources 
previously observed, for which the spectral index has not been determined. The empty circle 
marks the position of the 2.2 /xm peak (Lester et al. 1990). The newly detected features are 
labeled according to their positions (e.g. 39.92+55.9). 
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Fig. 3. — High angular resolution images of the H92q; hne emission (contours) throughout 
the velocity region where emission was observed in the central region of M82. The beam size 
is 0''6, shown in the bottom left corner of the first image. The contour levels are -3, 3, 4, 5, 

6 times the rms noise of ~0.12 mJy beam~^. The central heliocentric velocity is given above 
each image. The cross (+) indicates the position of the compact source 43.00+59.0, classified 
in this work as SNR (sec discussion in section 5), and the rotated crosses (X) indicate the 
HII regions; all of these arc compact sources observed in the continuum (McDonald et al. 
2002; Huang et al. 1994) and coincide with features observed in the H92a line. 
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Fig. 4. — Top) Contour images of the H92a line intensity throughout the velocity region 
where emission was observed at intermediate-resolution (0''9). The contour levels are —3, 
3, 4, 5, 6, 7, 8, 9 times the rms noise of ~0.07 mJy beam^^. The integrated line emission 
obtained from the low-resolution (2") line images is shown in gray scale. The gray scale 
image covers the range from to 6 Jy beam^^ m s^^. The marks indicate the position of the 
HII regions (crosses) and the SNRs (stars). The central heliocentric velocity is given above 
each image. The beam size is shown in the bottom left corner of the first image. Bottom) 
Integrated H92q; line emission at 0'.'9 angular resolution toward M82. Contour levels are 0.4, 
0.8, 1.2, 1.6, 2.0 mJy beam~^ km s~^. The cross indicates the 2.2 /xm peak (Lester et al. 
1990). The outer ellipse shows the extension of the inner 2.2 //m plateau. The ellipses shown 
correspond to models of Table 6 in row 1 and 5, respectively. 
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Fig. 5. — Top) Contour images of the H92« line at low angular resolution (2"). The contour 
levels are —3, 3, 5, 7, 9, 11, 13, 15, 17, 19 times the rms noise of ~0.07 mJy beam^^. The 
filled circle shows the position of G41.95+57.5 and the filled square the peak at 2.2 fim 
(Lester et al. 1990). The beam size is shown in the bottom left corner of the first image. 
Bottom) Integrated H92q; hne emission at 2" angular resolution toward M82. Contour levels 
are from 15.82 up to 158.2 by increment of 15.82 mJy beam~^ km s~^. The two line segments 
show the minor and major axis of a thin disk which would be observed with an inclination of 
81° and centered at the position of the 2.2 iim peak (Dietz et al. 1986). The cross indicates 
the 2.2 iim peak (Lester et al. 1990). 
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Fig. 6. — a) The integrated H92q! line emission (contours) superposed on the radio continuum 
at 8.3 GHz (gray-scale) at 2" angular resolution. Contour levels are drawn at 18 through 
148 mJy beam~^ km s~^ in steps of 26 mJy beam"^ km s~^. The gray scale for the continuum 
covers the range from 5.4 to 55 mJy beam~^. b) The integrated H53q; line emission (contours) 
superposed on the radio continuum image of M82 at 43 GHz (gray-scale) at 2" angular 
resolution. Contour levels are at 57, 114, 171, 228, 285 Jy beam~^ km s"^. The gray scale 
for the continuum covers the range from 2.3 to 34 mJy beam"^. c) Integrated intensity of 
H92q; hne (gray-scale) and integrated intensity image of the H53q; hne (contours). Contour 
levels are 57, 114, 171, 228 Jy beam"-^ km s~^. The gray scale covers the range from 8.4 to 
196 mJy beam^^ km s^^. The square marks the position of the 2.2 //m peak (Lester et al. 
1990) and the circle marks the position of the SNR G41. 95-1-57.5. 
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Fig. 7. — H92a and H53q; spectra toward the regions labeled in Figure 6c. The dashed curves 
are the observed intensities and the solid lines represent the fitted gaussians as a function of 
heliocentric velocity. 
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Fig. 8. — a) H92a velocity field at 2" angular resolution. Contour intervals are 16 km s~^, 
the iso-velocity contours > 200 km s~^ are in black, b) Spatial distribution for the velocity 
dispersion of the H92a profiles. The contours from 15 to 35 km are in black and above 
this level in white. The contour interval is 5 km s~^. The most elongated ellipse is a fit 
to a ridge with minimum velocity dispersion which corresponds to the transition between 
the inner and outer plateau in the 2.2 /im brightness distribution. In both images, the less 
elongated ellipse corresponds to the model of row 5 (Table 6). The two perpendicular line 
segments are the same as in Figure 5 (bottom). The cross indicates the center of the ring 
fitted by Achtermann and Lacy (1995). 
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Fig. 9. — Expected variation of radio continuum (left) and integrated line flux density (right) 
for four of the regions identified in Fig. 6c. The continuum fiux density has been computed 
using the model of multiple HII regions with a single component plus non-thermal emission 
arising from SNR. The integrated line emission curves correspond to the single- density models 
of Table 5. 
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Fig. 10. — H92a position-velocity diagram summed over the minor axis (angular resolution 
0''9), along: top) the main axis of the stellar bar with PA=70°. Contours are from 0.3 to 3.0 
in steps of 0.6 mJy beam^^, and bottom) the major axis of M82 determined on larger scales 
with PA=65°. Small circles marks the velocity of the stars as determined from [Ca II] hne 
(McKeith et al. 1993) and big circles show the velocity of the gas determined in Pa(lO) hne 
(McKeith et al. 1993). Contours are 0.1, 0.4 and 1.5 mJy beam-^ The beam is 2'.'5 x 0'.'9. 
The gray scales follow the contours. The slices were centered at Q;=09'*51'"43f6, 5=69°55'00". 
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Fig. 11. — Position-velocity image along the E ridge of broad profiles (see Fig. 8) at 150°. 
The angular resolution is 0''9. The offsets are relative to the major axis, NW to the left 
and SE to the right, the origin being at a =09'*51'"43f81, 5 =69°55'l."'l. Contours are -0.1, 
-0.05 mJy beam~^ (dashed) and from 0.05 to 0.5 mJy beam~^ in steps of 0.05 mJy beam~^. 



